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Abstract 


We  show  that  the  far  fields  generated  by  a  source  inside  or  near  a  micropar¬ 
ticle  can  be  obtained  readily  by  the  use  of  the  reciprocity  theorem  along 
with  the  internal  or  near  fields  generated  by  plane-wave  illumination.  The 
method  is  useful  for  solving  problems  for  which  the  scattered  fields  gener¬ 
ated  with  plane-wave  illumination  have  already  been  obtained.  We  illustrate 
the  method  for  the  case  of  a  homogeneous  sphere,  and  then  apply  it  to  the 
problem  of  emission  from  a  dipole  inside  a  sphere  near  a  plane  interface. 
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1.  Introduction 


Methods  for  obtaining  the  radiation  from  a  source  inside  or  near  a  micropar¬ 
ticle  are  needed  for  a  variety  of  applications:  e.g.,  in  modeling  the  fluorescence 
[1-4],  Raman  [5,6],  lasing  [7-9],  or  nonlinear  emission  [10,11]  from  molecules 
inside  or  near  scattering  objects  such  as  homogeneous  [12,13]  or  layered  [14] 
spheres,  spheroids  [15],  cylinders  [16],  microdisks  [8],  radially  inhomogeneous 
bodies  [17],  particles  with  complex  structures  [18],  or  particles  near  surfaces 
[19-33].  Methods  for  modeling  emission  from  polarization  sources  inside  par¬ 
ticles  are  also  needed  for  some  techniques  [34]  for  calculating  light  scattering 
by  inhomogeneities  inside  particles  [35,36]. 

The  problem  of  emission  from  a  molecule  or  other  point  source  inside  or 
near  a  microparticle  can  be  modeled  if  we  treat  each  point  source  as  a  time 
harmonic  (exp  {—iut)  time  variation)  dipole  p(rb)  at  r?,,  which  generates  an 
electric  field  E(ra)  at  Fa.  We  can  model  emission  from  molecules  with  nonzero 
emission  linewidths  by  integrating  the  electric  fields  over  emission  frequencies 
[3].  The  generated  fields  are  related  to  the  dipole  source  as  follows: 


Gii(ra,r(,)  Gi2(ra,rh)  Gi3(r„,rh) 

Pi(rb) 

Esiva)  . 

G'2l(ra5rb)  G22{^a,'f^b)  G23(^a,^b) 

_G3i(ra,r(,)  G32(ra,rb)  G33(ra,rb)  _ 

P2(rb) 

.  PsM  _ 

where  the  matrix  (typically  labeled  G(ra,r(,))  is  the  dyadic  Green  function 
[14,37],*  and  where  the  cu-dependence  of  E,  G,  and  p  is  suppressed.  We  have 

•Another  common  way  to  write  the  Green  function  relation  is 

E(r)  =  [  G(r,r')P(r')  dv', 

Jv 

where  P(r')  =  -zwJfr').  The  dipole  moment  p(r')  of  the  source  is  related  to  the  polar¬ 
ization  per  unit  volume,  P(r')  by 

p(r')  =  f  P(r')  dv'. 

Jv 

We  use  the  notation  of  individual  dipoles  because  we  have  been  modeling  radiation  from 
individual  molecules. 


1 


assumed  that  the  permeability  p  is  uniform.*  Although  equation  (1)  is  valid 
for  general  Tq  and  ti,,  in  this  report  we  treat  only  the  case  in  which  Tq  is  far 
from  the  particle,  and  is  inside,  on,  or  near  it. 

The  Green  function  of  equation  (1)  obeys  the  reciprocity  relation  [14,38,39],'^ 

G(r„,r6)  =  G^(r6,ra),  (2) 

where  T  indicates  transpose.  Given  a  solution  to  a  scattering  problem  with  a 
source  at  point  r^,  we  can  use  the  reciprocity  of  the  Green  function  to  obtain 
[13]  or  verify  [14]  Green-function  solutions  to  scattering/emission  problems 
for  sources  in  other  regions.  Beginning  with  the  Green  function  for  a  source 
near  a  sphere,  we  can  use  reciprocity  to  obtain  the  solutions  for  the  fields 
generated  by  an  incident  plane  wave  [37]. 

In  this  report  we  describe  a  simple,  reciprocity-based  method  for  obtaining 
the  far-held  Green  function  (the  Green  function  for  helds  far  from  an  object) 
for  a  source  inside  or  near  a  microparticle  or  other  scattering  object,  when 
the  solutions  for  the  helds  generated  by  an  incident  plane  wave  are  known. 
This  far-held  Green  function  differs  from  the  complete  Green  function  for  an 
emission  problem  in  that  several  elements  of  the  complete  Green  function 
are  not  specihed.  The  far-held  Green  function  is  often  all  that  is  required  be¬ 
cause  it  is  typically  all  that  can  be  detected,  for  example,  by  a  lens-detector 
system  located  far  from  the  particle.  A  main  beneht  of  the  approach  is  that 
the  far-held  Green  function  can  be  obtained  in  a  simple  manner  from  existing 
solutions  for  the  helds  generated  with  plane-wave  excitation.  Such  helds  have 
been  obtained  and  implemented  in  computer  codes  for  a  variety  of  scattering 
objects  (e.g.,  homogeneous  and  layered  spheres  and  spheroids,  spheres  with 
continuously  variable  refractive  index,  hnite  cylinders,  objects  with  axisym- 
metric  surfaces  described  by  Ghebyshev  polynomials,  cylinders,  and  particles 
on  or  near  plane  interfaces  [21-26]).  Another  beneht  of  the  approach  is  that 
we  can  transfer  some  of  the  understanding/intuition  developed  for  internal 
and  near  helds  of  spheres  [40-42]  and  cylinders  [42]  to  the  emission  problem. 
With  the  reciprocity  relations  for  plane  waves  and  far  helds,  the  intuition 
developed  with  ray-optics  or  other  methods  can  more  readily  be  used  to  as¬ 
sist  in  understanding  the  problem  of  emission  from  a  source  inside  a  particle 
[3].  We  can  also  use  known  solutions  for  plane- wave  incidence  to  validate 
the  far-held  limits  of  newly  developed  Green  functions  and  their  computer 
implementations . 

*The  assumption  of  a  uniform  permeability  is  valid  for  the  problems  we  want  to  model, 
which  are  at  optical  frequencies. 

^See  Chew  [14],  pp  410-411.  The  relation  for  regions  with  varying  is  G(ra,r(,)/i(rfc)  = 

G^(r6,ra)/x(ra). 


In  section  2,  we  show  how  the  Green  function  for  the  fields  far  from  a  mi¬ 
croparticle  can  be  obtained  from  expressions  for  the  fields  generated  inside 
or  near  a  microparticle  by  an  incident  plane  wave.  In  section  3,  we  illustrate 
how  this  method  works  for  the  case  of  a  sphere,  a  well-studied  particle  for 
which  the  Green  function  is  known.  In  section  4,  we  apply  the  method  to 
obtaining  the  Green  function  for  a  dipole  inside  a  sphere  on  or  near  a  con¬ 
ducting  surface,  a  problem  for  which  (so  far  as  we  know)  solutions  have  not 
yet  been  derived.  This  example  illustrates  how  readily  the  desired  expressions 
are  obtained  from  the  solution  to  the  plane-wave-incidence  problem.  Section 
5  summarizes  the  paper. 
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2.  Green  Function  from  Reciprocity 


We  assume  that  the  solution  to  the  problem  of  a  plane  wave  illuminating  the 
particle  is  available.  Our  goal  is  to  obtain  that  part  of  G(ra,r5)  required  for 
writing  the  fields  far  from  the  dipole  and  particle. 

We  begin  by  writing  the  plane-wave-illumination  problem  in  matrix  form; 


El{Tb) 

Fii(f6,F„)  Fi2(F6,Fa)  Fi3(rfc,Fa) 

EoBo  ■  il(Fa) 

E2{n) 

= 

^2i(r6,ra)  E22(ri„ra)  E23(rb,ra) 

EqBq  •  l2(Fa) 

Esin)  _ 

.  F3i(Fb,Fa)  E32(rb,ra)  F33(r6,Fa)_ 

_  EoBo  ■  i3(ra)  . 

where  Ei{rt),  E2{rb),  and  Esin)  are  the  field  components  at  position  n  inside 
or  near  the  particle,  and  the  incident  plane-wave  field  is  given  by 

E'""(r)  =  EoeoE^",  (4) 


where  k  is  the  propagation  vector  and  the  unit  vector  Cq  is  perpendicular  to 
k.  We  write  the  unit  vectors  as  ij(ra)  to  emphasize  that  each  is  evaluated 
at  Tfl.  An  example  of  Eij  for  axisymmetric  particles  analyzed  in  a  spherical 
coordinate  system  is  given  in  appendix  A. 


The  incident  plane  wave  can  be  generated  to  any  degree  of  accuracy  by  a  time- 
harmonic  dipole  polarization  source,  p(ra)  =  eop(ra),  with  60  perpendicular 
to  Fa,  when  the  dipole  is  sufficiently  far  from  the  object  (i.e.,  |ra|  >  A, 
and  I  Fa  I  is  many  times  larger  than  the  object).  (The  surface  on  which  the 
outgoing  boundary  conditions  are  applied  for  the  Green-function  solutions  is 
many  times  further  from  the  object  than  1fq|.)  Near  the  particle,  the  field  of 


the  dipole  is  [43] 


E(f) 


47r  |F-Fa| 


p(ra)eo. 


(5) 


Comparing  equations  (4)  and  (5),  we  see  that  the  amplitude  of  the  plane 
wave  generated  by  the  source  p(Fa)  is  approximately 


up' UP 

^  ^ 
47rr„ 


(6) 


where  |f  —  Fo|  in  the  denominator  is  replaced  by  fa  =  | Fa |  because  the  points 
F  are  near  the  particle  and  because  |f|  <?C  [fo]. 
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Using  equation  (6)  for  Eo,  we  can  rewrite  equation  (3)  as 

jUi(rfc)  Gii(r(,,ra)  Gi2(r6,ra)  Gi3(r6,ra)  Pi(ra) 

E2{n)  G2i(rb,ra)  G22(r6,ra)  G23(r6,r„)  P2{ra)  >  (7) 

_  Esirb)  _  _G3i{rb,Ta)  G32(rb,ra)  G33(rb,ra)  _  ^3(0. 

where  ^ 

Gij(rb,ra)  =  ^^;^Eij(rb,ra).  (8) 

Equations  (7)  and  (8)  provide  the  Green  function  for  the  fields  inside  an 
object  excited  by  a  dipole  far  from  the  object.  Then,  usmg  the  reciprocity 
relation,  equation  (2),  we  obtain  the  desired  elements  of  G(ra,rb)  generated 
by  a  source  at  inside  or  near  the  particle: 

Ei(r„)  1  Gii(rb,r„)  G2i(r6,ra)  G3i(r6,r„)  pi(rt) 

E2(ra)  Gi2(r6,ra)  G22(rb,ra)  G32(rb,ra)  P2{n)  •  (9) 

_  E3{Ta)  _  _  Gi3(rb,ra)  G23(r6,ra)  G33(rb,ra)  _  _  P3(ri.)  _ 

It  must  be  emphasized  that  the  above  Green  function  is  valid  only  for  field 
points  Fa  far  from  the  particle.  For  Fq  close  to  the  particle,  some  or  all  of  the 
Gji  would  have  additional  terms  that  would  decay  with  distance  from  the 
particle. 

If  a  spherical  coordinate  system  is  chosen,  then  f^  is  far  from  the  particle, 
on  a  line  going  from  the  origin  in  the  — k  (or  U  direction),  and  pi(ra)  for  the 
incident  plane  wave  is  zero;  thus,  the  plane-wave-illumination  problem  can 
be  written  as 

Er{rb)  U[i{rb,Ta)  Fi2(f6,F„)  Fi3(Fb,Fa)  0 

Eeirb)  =  U2i{Tb,ra)  F22{rb,ra)  F23{rb,Ta)  EoBo-ie  >  (^0) 

Etbipb)  ,  G^(Fb,Fa)  E32(Fb,  Fq)  E33(Fb,Fa)  _  _  EqGq  •  i(j)  _ 

where  we  write  the  Fii{rb,Ta)  as  U[j{rb,ra)  to  emphasize  that  they  are  un¬ 
specified  and  unneeded  for  the  internal/near  fields  at  r^.  The  Fi2(r6,Fa)  and 
Fj3(Fb,  Fa)  are  known.  The  desired  Green  function,  obtained  with  equations 
(8)  and  (2),  is 

0  1  r  U^i{rb,Ta)  Ui{Tb,ra)  Ug{rb,ra)  Prin) 

Ee{ra)  Gi2(rb,Fa)  G22(rb,ra)  G32(rb,Fa)  Pein)  ’ 

F4>{T^a)  .  .  Gi3(Fb,Fa)  G23(rb,  Ta)  G33(Fb,  Fa)  _  _  P(^(rb)  _ 
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where  Gii(rfc,r«),  G2i{vb,r^),  and  G3i(r5,ra)  are  written  as  f/jl(rfc,r„),  simi¬ 
lar  to  the  Ulin,  Ta)  of  equation  (10),  because  the  values  are  unspecified  and 
not  needed  for  the  far-field  solutions  at  Fq.  The  above  expression  is  valid  for 
complex  u,  and  so  may  be  useful  for  treating  problems  in  terms  of  quasinor¬ 
mal  modes  [44]. 

The  relations  given  in  equations  (9)  and  (11)  are  key  results  of  this  report. 
In  particular,  in  the  spherical  coordinate  system  commonly  used  for  parti¬ 
cle  scattering  problems,  equation  (11)  indicates  how  the  known  internal  field 
solutions  (the  i^nm)  specify  the  far  fields  from  an  arbitrary  dipole  source 
inside  the  microparticle.  The  understanding  of  the  spatial  variations  of  the 
Fnm  (proportional  to  the  Gmn)  for  some  microparticles  (spheres  [40-42]  and 
cylinders  [42])  can  now  more  readily  be  used  to  visualize  and  understand 
the  emission  problem.  Distributions  of  fluorescence  collected  from  oriented 
dipoles  inside  a  sphere  have  been  previously  shown  [3].  Although  the  previ¬ 
ous  study  noted  that  the  emission  problem  was  related  to  the  incident  field 
problem  by  reciprocity,  the  specific  relations  shown  in  equation  (11)  were  not 
known  at  that  time. 


3.  Example:  Fields  from  a  Source  in  a  Homogeneous 
Sphere 


To  illustrate  the  above  method  for  a  well-studied  particle  (a  homogeneous 
sphere),  we  compare  the  fields  at  Fa  on  the  -z  axis  far  from  the  particle 
generated  by  a  dipole  inside  the  sphere  and  lying  in  the  cj)  =  0  (i.e.,  x  —  z) 
plane.  In  this  case,  we  use  the  spherical  coordinate  system;  i.e.,  ii  =  ir  is 
the  unit  vector  in  the  radial  direction,  12  =  ie,  and  is  =  i^.  Expressions 
for  the  fields  emitted  by  a  source  inside  a  sphere  are  well  known  [12].  The 
fluorescence  collected  from  a  dipole  inside  a  sphere,  calculated  as  a  function 
of  dipole  position,  has  been  illustrated  elsewhere  [3]. 

3.1  Using  Reciprocity 

For  points  in  the  0  =  0  plane,  Fis{rb,ra),  E23(rb,ra),  and  F32(r(,,ra)  are  zero 
(see  app  B),  and  equation  (3)  reduces  to 

Er{rb)  U[^{rb,Ta)  Fi2{Tb,ra)  0  0 

Ee{rb)  =  Ui[{Tb,ra)  F22{rb,ra)  0  EoBo-ieiT^a)  ■  (^2) 

_  E^in)  _  [  F3^i(rb,r„)  0  F33(r6,  r„)  J  [  •  i^(r„)  _ 

The  desired  Green  function,  as  in  equations  (9)  and  (11),  is  of  the  form 

Er{Xa)  F(^(r6,rtt)  f/;^(rb,ra)  U^(rb,ra)  Pr(^b) 

Ee{ra)  G'i2(r6,ra)  G'22(r6,r„)  0  ^^(ri,)  •  (13) 

E^{Ta)  _  _  0  0  G33(rf„ra)  _  _P4,{rb)  _ 

We  consider  the  case  of  a  0-polarized  dipole  in  which  Pi  =  P2  =  0,  and 

p^{rb)  =  P4,-  Using  the  Green  function  from  reciprocity  (eq  (13))  and  equation 
(8),  we  obtain  Ee{Ta)  —  0,  and 

gikra 

E^iVa)  =  ra)P^  ^  fa),  (14) 

where  k  ■  Fq  =  —kva  because  f^  is  on  the  —2:  axis. 

Using  the  F33  given  by  equation  (A-6)  with  m  =  1,  we  obtain  from  equation 
(14) 
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£«(r,)  = 

-T. 


AkVa 


A-KTa 


(15) 


[rtk»n)] 


where  from  equations  (4.3)  and  (4.5)  of  Barber  and  Hill  [15] 
..  2n+l  ^ 


Celn  ^ 

rl  —  _7"+l 

Uoln  ^ 


n{n  +  1)  xjn{r]x)[xh^n\x)y  -  [T]xjnivx)]'xhly\x) 

2n  +  l  rji 


n{n  +  1)  T]‘^xjn{riA[^h^rl\x)]'  -  [7]xjn{vx)]'xh\t’{x 


,(1)/ 


(17) 


where  x  is  the  size  parameter  of  the  sphere  and  r]  is  its  refractive  index.  With 
equations  (16)  and  (17),  equation  (15)  can  be  written 


_  ^  (2n  +  1) 

”  kra  n  47r  n{n  +  1) 

z”  d 


+ 


smfe 

r]^xjnivx)[xhyi\x)]'  -  [pxjnivx)]' xh^n\x) 


*”'^^jn(hfc?’6)^-Pn(cOS  Ob) 
XJn{Vx)[xh^A  {x)]'  -  [T]Xjnivx)]' Xh^n\x) 

P^(C0S  Ob) 


(18) 


3.2  Using  the  Complete  Green-Function  Solution 

We  verify  our  results  by  comparing  them  with  results  obtained  in  the  tra¬ 
ditional  boundary- value  fashion:  i.e.,  starting  out  with  a  radiating  dipole 
within  the  sphere,  satisfying  the  boundary  conditions  at  the  sphere  surface 
and  finding  the  resulting  scattered  far  field  in  the  =  tt  direction.  Using 
the  complete  Green-function  solution  for  a  source  inside  a  sphere  (see  app  B) 
with  P3(rfc)  =  and  pi  —  P2  =  0,  we  obtain  the  scattered  field  coefficients 
as 


fv^  =  -  App^ — 


^[^^{pkrb) 


jn{rjA[^'^n\x)\  -  [pxin{'nx)]h!'n\xy 


by  using  equations  (B-4)  and  (B-7),  and 


(19) 
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(20) 


k 


rji^  ■  Nu^{T]kTb) 


^  '^^^7rxrj2j^(^r]x)[xh^n\x)]'  -  [rjxjnivx)]' h^n\xy 

by  using  equations  (B-5)  and  (B-8). 

Prom  the  definitions  of  M.i'^{r]kTb)  and  Nu^{r]krb),  the  and  gu^  simplify 
to  fZin  =  9?mn  =  0>  ^nd 


/, 


Lo^gikp^ 


jn{vkrb)^Pn  (cos  Ob) 


emn 


Q 

9omn 


TTX 


jn{vx)[xh^n\x)]'  -  [r]Xjn{gx)]'hy\xy 


-Lo'^gp^, 


±,ri  (vkn)]SliSB^ 


irxTb  r)^j„(rix)\xhn\x)]'  -  {qxjnirjx)]' h^\x) 


(21) 


(22) 


Now,  by  restricting  the  field  points  Va  to  the  z  axis  where  the  contributions 
for  m  1  are  zero,  and  restricting  the  source  points  rt  to  the  (^  =  0  plane, 
noting  that 


-Pi  (cos  0a) 


sin  6  a 


Oa~TT 


0a=7r 


=  (-1) 


(-1) 


n+Mn±l) 


,  n 


(’»+ 1) 


(23) 

(24) 


(because  0^  =  tt  at  the  field  point),  and  using  the  far-field  expansion  of  the 
spherical  Hankel  function. 


hy\kr)  = 


(n+i) 

kr 


Akr 


(25) 


we  obtain 


Mei„(*:r)|,,,  =  i 


kr 


i(j>i 


N  (fcrll  - 

Then  using  equation  (B-6),  we  obtain  Eg  =  0,  and 


(26) 

(27) 


6*'=’'“  ^  n{n  +1) 


In 


q(«+l) 


feln  “t  ^  9oln 


(28) 
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which  is  equivalent  to 


gifcra  (2n  +  1) 


i("+i)j„(r7fcrfc)— P„Hcos  9b) 


An  n{n+l)  xjn{vx)[xhK\x)]' -  [r]xjn{vx)]' xh\^\x) 


+  krbdn^'^^^^  sinOb 

TJ-^XJnMlxhi^Hx)]'  -  [vxjn{vx)]'xhi^^\x) 


The  expressions  for  the  electric  field  that  we  obtained  using  reciprocity  (eq 
(18))  and  using  the  complete  solution  (eq  (29))  are  the  same. 


4.  Example:  Fields  from  Source  in  Homogeneous 
Sphere  near  a  Plane  Conducting  Interface 


Fluorescence  and  Raman  emission  have  been  used  in  characterizing  particles 
on,  inside,  or  near  a  surface,  e.g.,  a  biological  cell  or  spore  on  a  filter,  or  a 
contaminant  particle  on  a  silicon  wafer.  Solutions  have  been  described  for 
the  fields  scattered  by  a  sphere  in  close  proximity  to  a  plane  interface  and 
illuminated  with  a  plane  wave  [21-26].  However,  as  far  as  we  know,  solutions 
for  the  fields  emitted  from  a  dipole  inside  a  sphere  on  a  plane  surface  have  not 
been  described.  Here  we  use  the  reciprocity  theorem  and  a  known  solution 
for  the  fields  generated  in  a  sphere  with  plane-wave  excitation  to  write  the 
solution  for  the  far  fields  generated  by  a  dipole  emitting  inside  a  sphere  on  or 
near  a  plane  surface.  To  obtain  the  required  for  equation  (3),  we  use  the 
derivation  of  Videen  [26]  who  presents  a  solution  to  the  fields  of  a  particle  on 
or  near  a  perfectly  conducting  plane  surface.  The  vector  spherical  harmonics 
used  in  that  article  [26]  and  in  this  section  are  normalized: 


-  j,zi''Hkr)^Pr(cose)e‘’"*, 


sin  9 


f^zlf\kr)n{n+l)P;ri^os9)e^^‘^  +  9^-  [rzi^\kr)]  ±P^{cos9)e 


dO 
kr  dr 


(30) 


dO' 


im<j) 


Although  it  is  usually  not  appropriate  to  use  multiple  definitions  of  the  vector 
spherical  harmonics  in  one  paper,  here  the  multiple  definitions  help  illustrate 
the  generality  of  the  approach  stated  in  section  2. 

We  solved  the  plane-wave-incidence  problem  by  expanding  the  incident  plane 
wave,  the  scattered  fields,  and  fields  interior  to  the  sphere  in  vector  spher¬ 
ical  harmonics.  In  addition  to  these  fields,  there  is  an  interaction  field  that 
scatters  from  the  sphere,  reflects  from  the  plane  surface,  and  illuminates  the 
sphere  again.  We  determine  the  field  coefficients  by  forcing  the  boundary 
conditions  at  the  interfaces  of  the  sphere  and  plane  surface  to  be  satisfied 
simultaneously. 

The  fields  inside  the  homogeneous  sphere  are  expanded  as  [26] 
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(31) 


Efivkr)  =  Y,  (,,fcr), 

n,m 

where  &re  the  interior  field  coefficients.  These  coefficients  are  expressed 
in  terms  of  the  known  scattering,  interaction,  and  incident  field  coefficients, 
blil,  and  as 

(^nijnivkr)  =  a^nljn  {^a)  +  b^hn  {ka)  +  {ka) ,  (32) 

^(^nmJnivkr)  =  Clnmjn  {ka)  +  {ka)  +  C^n^jn  {ka)  .  (33) 


We  can  express  equations  (31)  to  (33)  in  the  Green-function  formalism  of 
section  2  by  writing  out  the  vector  spherical  harmonics  in  equation  (31),  and 
writing  the  fields  in  terms  of  the  Fij  of  equation  (3)  as 


i^i2(r6,r„)  = 


Fi3{rb,ra) 


T22(r6,r„) 


T23(r6,ra)  = 


F32{rb,ra)  = 


F33{rb,ra)  = 


n,m  I  0 

E  +  l)P”(cos  (35) 

n,m  TjICTf) 

n,m  Sm  &b 

E  e^r-^^llHnkn)PT(cos  ^)e— *  +  (37) 

sm  &b 

E  -4)f*'4"(-?*:>'^)  J-fr(cos9.)e‘--  +  (38) 

n,m  ^ 

b4‘>(.'cn)]  £fPMy-\ 

E  -e(«™4'>(»*:n)^Pr(cos»i,)e‘”*‘  +  (39) 

n,m  ^ 

rjkvbdrb  L  "  ^  'i  smdb 


where  the  superscripts  TE  and  TM,  on  the  internal  field  coefficients,  refer 
to  the  polarization  state  of  the  incident  plane  wave.  The  desired  Green  func¬ 
tion  is  then  given  directly  by  equation  (11)  with  Gij{Tb,Ya)  obtained  from 
Fij(r6,ra)  as  described  in  equation  (8).  Thus,  the  scattered  field  components 
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resulting  from  a  dipole  within  a  sphere  near  a  perfectly  conducting  surface 
are  obtained  with  equation  (31). 

Although  the  above  was  derived  for  the  particular  problem  of  a  sphere  above 
a  surface,  the  preceding  expressions  for  the  Fij{rb,ra)  and  the  Gtj(rh,ra)  are 
valid  for  any  particle  for  which  the  internal  field  coefficients  of  equation  (31), 
elli  and  are  known.  For  a  dipole  outside  but  near  the  particle,  equations 
(34)  to  (39)  may  also  be  used,  where  the  internal  field  coefficients  *  are 
now  replaced  by  the  scattered  field  coefficients,  and  the  Bessel  functions  are 
replaced  by  the  appropriate  Hankel  functions.  The  particular  case  of  a  sphere 
near  a  substrate  is  slightly  more  complicated  when  the  dipole  is  outside  the 
sphere,  because  the  interaction  field  must  also  be  included  with  the  scattered 
field.  Finally,  although  we  have  used  normalized  vector  spherical  harmonics, 
these  equations  may  also  be  applied  to  coefficients  derived  with  unnormal¬ 
ized  vector  spherical  harmonics  by  a  simple  replacement  of  the  normalized 
associated  Legendre  polynomials,  P^{cos6b),  with  the  associated  Legendre 
polynomials,  P^{cos6b). 
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5.  Summary 


This  report  is  based  on  the  following  three  observations:  (1)  The  Green  func¬ 
tion  relating  a  source  and  a  scattered  field  obeys  a  reciprocity  relation  [14,38]. 
(2)  Solutions  for  the  fields  inside  or  near  a  variety  of  scattering  objects  have 
been  obtained  for  plane- wave  excitation  but  not  more  general  sources.  (3)  For 
many  problems  in  which  a  molecule  or  polarization  source  emits  inside  or  near 
a  scattering  object,  it  is  only  the  far  fields  that  are  measured  or  are  of  in¬ 
terest;  for  these  problems,  a  far-field  Green  function  that  only  specifies  the 
fields  far  from  the  particle  is  sufficient. 

The  key  point  of  this  report  is  that  reciprocity  and  known  solutions  for  the 
fields  generated  in  a  particle  by  plane-wave  incident  fields  can  be  readily  used 
to  find  the  far  fields  emitted  by  a  source  inside  or  near  the  particle.  Although 
only  a  partial  Green  function  is  known  from  the  plane-wave-incidence  prob¬ 
lem,  that  partial  Green  function  and  reciprocity  are  sufficient  for  specifying 
the  far  fields.  These  key  results  are  given  in  equations  (9)  and  (11). 

In  section  4,  we  illustrate  the  technique  by  applying  it  to  a  particle  for  which 
the  solution  is  well  known:  a  homogeneous  sphere  in  a  homogeneous  medium. 
We  then  demonstrate  the  power  of  this  technique  by  applying  it  to  a  more 
complicated  problem  for  which  the  fields  emitted  by  a  dipole  in  the  parti¬ 
cle  are  not  known:  a  dipole  located  within  a  sphere  near  a  plane  interface. 
Although  the  solution  was  derived  for  this  particular  system,  the  equations 
given  in  section  4  can  be  applied  to  any  system  for  which  the  internal  or  scat¬ 
tering  coefficients  of  an  expansion  in  vector  spherical  harmonics  have  been 
derived. 

Another  benefit  of  understanding  these  reciprocity  relations  of  the  inter¬ 
nal  fields  generated  by  incident  plane  waves  is  that  the  understanding  of 
the  internal  intensity  patterns  of  particles  (obtained  from  ray-optic  analyses, 
comparisons  with  Fabry-Perot  cavities,  etc)  can  be  used  to  help  develop  an 
understanding  of  emission  patterns  from  sources  inside  the  particle  [3]. 
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Appendix  A.  Example:  for  Axisymmetric 

Particles 


If  the  scattering  solution  is  found  for  an  axisymmetric  particle  with  a  spher¬ 
ical  coordinate  system  and  spherical  wave  functions,  and  r;,  is  in  the  0  =  0 
plane,  then  F12,  F22,  ^32  are  given  by  the  summation  over  n'  of  Ef\ 
E'^^  of  equation  (3.10)  of  Barber  and  Hill  [15],  i.e., 


Fi2 

E22 

F32 


V  n  (n  +  1)  cos  sin  Obdemn, 

^  rjkn  smOb 


dn  ivkrb)  cos  m06 


mP™  (cos  6b) 


+  - 3 —  [^"6in(^^^6)]  cos  7720  (cOS  0;,)  demm 

Tjkrb  dvb  u" 

Y  ivkrb)  sinm0fc^P,("(cos  6b)Comn 


1  d 


mP^  (cosdb) 


hjni^krb)]  sinm0fc - 

rjkrb  dvb  sm  Ob 


(A-1) 

(A-2) 

(A-3) 


where  the  Cemn,  Comn,  demn,  and  damn  are  the  internal  field  coefficients  of  the 
particle.  The  Fn  are  not  specified  in  the  plane-wave-incidence  problem.  The 
Pi3,  -^23,  -^33  are  given  by  the  summation  over  n'  of  P‘"^  E'^\  E'^^  of  equation 
(3.11)  of  Barber  and  Hill  [15]: 


Pi3  — 
F23  = 

F33  = 


V-  /  ,  1  ^  F  ivkrb)  ^  Pn  (cos  6b)  ^  ^ 

/  Ti  {n  4*  1) - z  sin  UKpi,  .  ^  sin  ui^domn-) 

4;  r}krb  sin  06 


Y-jn  {rjkrb)smm(f)b 


mP^  (cos  6b) 


+  Injnivkrb)]  sinm06^P"(cos  6b)d  omni 

Y  ~dn  {vkn)  cos  m4)b-^Pn  {cos  6b)Cemn 

n^m 

+  [rbjniilkrb)]  cosm(t)b^^'' damn- 

-qkrbdn  sm06 


(A-4) 

(A-5) 

(A-6) 
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If  the  particle  has  spherical  symmetry,  these  equations  can  be  simplified, 
since  only  the  m  =  1  terms  contribute  to  the  scattered  fields.  The  reduced 
equations  for  Fi2i  F22,  T32  are  given  by  the  summation  over  n  of  ,  Eg  , 
E'^^  of  equation  (4.33)  [15],  and  F13,  jp23i  F33  are  given  by  the  summation 
over  n  of  El'^\  Ef\  E'f  of  equation  (4.34). 


Appendix  B.  Fields  from  a  Dipole  Inside  a  Sphere 


The  Green  function  for  a  source  inside  a  sphere  has  been  described  [12  14,37]. 
The  total  internal  field  generated  by  the  source  inside  the  sphere  is  [3] 

E^{rikr)  =  E^{r]kr)  +  E'‘^(7?fcr).  (B-1) 

Here,  E"  {rjkr)  is  the  electric  field  at  r  radiated  by  arbitrary  polarization 
sources  at  r'  (where  r  >  r')  in  a  homogeneous  region  of  refractive  index  rj, 

OO 

E^{rikr)  =  ^  Di'[cu^M.i'^{r]kT)  +  di>^'Ni'^{r]kr)].  (B-2) 

W=1 

The  superscripts  on  the  vector  spherical  harmonics  refer  to  the  kind  of  radial 
function  (first  or  third)  used  in  the  expansion,  and  v  represents  the  spherical 
harmonic  triple  index  a,  m,  n,  where  a  is  even  or  odd,  n  is  the  mode  number, 
and  m  is  the  azimuthal  mode  number.  The  normalization  constant  is 

n  —  ^)(^  ~  ^)'  fB-3) 

4n(n  +  l)(n  +  m)!  ’ 

where  is  equal  to  1  for  m  =  0  and  equal  to  2  for  m  >  0. 

The  field  expansion  coefficients  cu^  and  du^  are  determined  by  the  strength, 
position,  and  orientation  of  the  source  polarization  p(rfc)  as 

p(rf,)  •  Mi/^(77fcrb),  (B-4) 

TT 

du^  =  p(rj  ■  Ni/^(7]/crt),  (B-5) 

TT 

where  the  superscript  H  indicates  a  homogeneous  region. 

The  induced  electric  fields  outside  the  sphere  are  expanded  as 

OO 

E®‘^(A;r)  =  ^  Du[fi>*^M.u^{kr)  +  gu^'Nu^{kr)],  (B-6) 

i/=i 

where  and  gu^  are  termed  the  “scattered”  field  expansion  coefficients 
because  of  their  similarity  in  equation  (B-6)  to  the  scattering  coefficients  of 
the  usual  scattering  problem.  The  G  in  the  superscripts  of  the 
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and  gu^  differentiates  these  internal  and  scattered  coefficients  from  the 
field  coefficients  used  with  other  incident  fields. 

The  scattered  field  coefficients  are 


rjxjn{vx)[xh^n\x)]'  -  r][r]xjn{vx)]' xh^n\x) 

i 

rfxjn{'nx)[xh!'n\x)]  -  (x) ’ 


(B-7) 

(B-8) 


where  x  is  the  size  parameter  and  g  is  the  refractive  index  of  the  host  sphere, 
{gx)  is  the  spherical  Bessel  function,  and  (x)  is  the  spherical  Hankel 
function  of  the  first  kind.  All  derivatives  (denoted  by  the  primes)  are  with 
respect  to  the  argument. 

In  our  earlier  work  describing  the  fields  from  a  dipole  inside  a  sphere,  the 
sign  of  the  coefficients  of  the  induced  internal  transverse  magnetic  fields  (the 
in  eq  (A9)  of  Hill  et  al.  [3]  and  in  eq  (12)  of  Hill  et  al.  [34])  should 
have  been  negative.  These  coefficients  are  not  used  here,  nor  were  they 
used  in  the  previous  references  [3,34].  However,  the  sign  error  is  important 
for  anyone  verifying  the  above  expressions. 
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